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Abstract 

i . 
O ■ 

O ■ .D-dimensional cosmological model describing the evolution of a perfect fluid with neg- 

ative pressure (x-fluid) and a fluid possessing both shear and bulk viscosity in n Ricci-flat 
spaces is investigated. The second equations of state are chosen in some special form of 
metric dependence of the shear and bulk viscosity coefficients. The equations of motion 
are integrated and the dynamical properties of the exact solutions are studied. It is shown 
■ the possibility to resolve the cosmic coincidence problem when the x-fluid plays role of 
quintessence and the viscous fluid is used as cold dark matter. 

(N ■ 
> ■ 

O ! 1 Introduction 

o ■ 

The most probable candidate for the so called quintessence matter responsible for the current 
phase of the accelerated expansion of the Universe is a A term o more generally an exotic 
x-fluid - perfect fluid with negative pressure satisfying a linear barotropic equation of state 
(see, for instance, |] and refs. therein). To describe the present stage of evolution this x-fluid 
is to be added to the normal matter, which mainly consist of cold dark matter. For instance, 
ACDM model describes the flat Friedman- Robertson- Walker (FRW) Universe filled 

with a mixture of quintessence represented by the cosmological constant A and cold dark 
^ . matter in the form of pressureless perfect fluid (dust). 

However, the flat FRW cosmologies with a x-fluid and a normal perfect fluid are not free 
from some difficulties. One of them is the so called coincidence problem ||]: to explain why 
the quintessence density and the normal matter density are comparable today, one has to 
tune their initial ratio very carefully. The problem may be ameliorated by replacement the 
x-fluid by the so called quintessence scalar field - homogeneous scalar field Q slowly rolling 
down with some potential V(Q) §,[§]. For instance, the potential V(Q) = M i+a Q~ a , a > 0, 
leads to the so called "tracker" solution, which "attracts" solutions to the equations of motion 
before the present stage for very wide range of initial conditions. Another way for resolving 
the problem was proposed in [[?] . The idea is to use a fluid with bulk viscosity in combination 
with a quintessence scalar field. 

The aim of this paper is to show that the cosmic coincidence problem may be resolved 
by using a x-fluid as the quintessence and a viscous fluid as the normal matter. In section 2 
we describe the general model and get basic equations. To integrate the equations of motion 
we choose the so called " second equations of state" , which provide us with the dependence 
of the shear and bulk viscosity coefficients on time, in some special form of their metric 
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dependence. The exact solutions for both isotropic and anisotropic case are obtained in 
section 3, where their dynamical properties are discussed. 



2 The general model 

We assume the following metric 

n 

ds 2 = -e 2 ^ (t) dt 2 + £ exp[2^(t)]d^, (2.1) 
i=i 

on the .D-dimensional space-time manifold 

M = R x Mi x . . . x M n , (2.2) 

where dsf is the metric of the Ricci-flat factor space M; of dimension d h ^(t) and x l (t) are 
scalar functions of the cosmic time t. a, L = exp[x l ] is the scale factor of the space Mj and the 
function j(t) determines a time gauge. The synchronous time t s is defined by the equation 
dt s = exp[ 7 (*)]dt. 

We consider a source of gravitational field in the form of 2-component cosmic fluid. The 
first one is a perfect fluid with a density p^(t) and a pressure p^(t). The second component 
supposed to be a viscous fluid. It is characterized by a density p^ 2 \t), a pressure p^(t), 
a bulk viscosity coefficient ((t) and a shear viscosity coefficient i](t). The overall energy- 
momentum tensor of the cosmic fluid reads 

7? = (p (1) + p (2) ) u»u v + (p {1) + p {2) - CO) P? - 2?7<, (2.3) 

where w M is the D-dimensional velocity of the fluid, 9 = denotes the scalar expansion, 
Py = &v + u ^u u is the projector on the (D — 1) -dimensional space orthogonal to u^, = 
\ (u a -[3 + U/3-a) P afl pP — (D — l^OPjf is the traceless shear tensor and p, v — 0, 1, . . . , D — 1. 

By assuming the comoving observer condition = ^Qe -7 ^, the overall energy-momentum 
tensor may be written as 

T M = T Mi) +T M2) ) (2.4) 

where 

(T^) =diag(-p(V 1 ),...,p( 1 )), (2.5) 
(T^) = diag(-,ro M 2) #, • • • ,pf } <t), (2-6) 

ki,k = l,...,di for i = l,...,n. Here pf 1 denotes the effective pressure including the 
dissipative contribution of the viscous fluid in the factor space described by the manifold 
Mj. It reads 



p f > = p (2) _ e -7 



C7o + 2 V (i 



7o 

x 



V £>- 1 



(2.7) 



where 

n 

l0 = J2d lX \ (2.8) 



=i 
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Furthermore, we assume that the barotropic equations of state hold 

p(«) = (l - h {a) ) p (a \ a = 1,2, 



(2.9) 



where the are constants such that h^ 1 ' ^ hP^. 

The Einstein equations R£ — \5„R = ^Tjf, where n 2 is the gravitational constant, can 
be written as R$ = k 2 (T^ - [T/(D - 2)]<J£). Further, the equations R° - \5° Q R = k 2 T ° and 
Rl = k 2 (T^ — [T/(D — 2)]5%), where a,b = 1, . . . , D, give the following equations of motion 



7o ~^2di(i 



i\2 



2 ^ e 27 (p(D +p (2) 



i=l 



& + (70 - i)^ 



D-2 



a=l 



D-2 



A fo-2r]( x 



7o 



D - 1 



(2.10) 
, (2.11) 



i = I, ... ,n. 



The energy conservation law S/^Tq = for a viscous fluid described by a tensor given 
by equation ( |2.6| ) reads 

n 

(2.12) 



(2) 



+ J2d i x i (p {2) +P? ] )=0. 



i=l 



Owing to the constraint S/^Tq = for the overall energy-momentum tensor given by equation 
(|2.3| ) the similar energy conservation law is valid for the perfect fluid 



p(i) + ( p (i) +p (i))£^ = . 



(2.13) 



i=l 



Taking into account equation 



ones easily integrates equation (|2.13| ). The result is 
pW = Ae^ hW - 2 ^, (2.14) 



where A is a positive constant. We note that the contribution of the perfect fluid component 
with = 2 to the overall energy momentum tensor is equivalent to the presence of A-term 
with A = k 2 A. 

Furthermore, by using equations Q2.10|) and (|2.14j ) the presence of the densities and 
p( 2 ) in equations ( [2.11 ) can be cancelled. Thus, we obtain the main governing set of equations 



& + (70 - i)^ 



k 2 A 



hi 1 ) - hP) , m , J,(2) 
e [/i( 1 )-2] 70 +2 7 , " 



D-2 



2 7 

k e' 



D-2 



7o - 2rj [x l - 



2{D-2) 

7o 



7o -E^(^ 



i\2 



i=l 



D - 1 



(2.15) 



We use an integration procedure which is based on the n-dimensional Minkowsky-like 
geometry. Let IR n be the real vector space and e 1; . . . , e n be the canonical basis in IR™ (i.e. 
ei = (1, 0, . . . , 0) etc). Let us define a symmetrical bilinear form (, ) on IR™ by 



(2.16) 
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Such a form is non-degenerate and has the pseudo-Euclidean signature (— , +,...,+) || . 
With this in mind, a vector y £ IR n is timelike, spacelike or isotropic respectively, if (y, y) 
takes negative, positive or null values respectively and two vectors y and z are orthogonal if 
(y, z) = 0. Hereafter, we use the following vectors 



u 



x 1 {t)ei + . . . + x n {t)e n , 
u 1 e 1 + . . . + w n e n , u l = 



D — 2 



(2.17) 
(2.18) 



where the covariant coordinates Uj of the vector u are introduced by the usual way. Moreover, 
we obtain 

D - 1 



(*> *) = X) d i{ ±% f ~ 7o> ( u > x )=7o, (u,u) 



i=i 



D-2 



;2.i9) 



Thus, using equations (|2.17 )-( p.l9 ) we rewrite the main governing set of equations in the 
following vector form 



x + (7o - 7)x 



k 2 A (h^ - h^) e^ C1) - 2 ^+^ + — (x,x) 



u 



- «v 



(2.20) 



Moreover, the density can be expressed as 



P 



(2) 



(x, x) 27 _ (1) 
2^2 P • 



(2.21) 



Now we summarize thermodynamics principles. The first law of thermodynamics applied 
to the viscous fluid reads 

TdS = d (p (2) V) + p (2) dV, (2.22) 



where V stands for a fluid volume in the whole space Mi x . . . x M n , S is an entropy in the 
volume V and T is a temperature of the viscous fluid. By assuming that the baryon particle 
number Nb in the volume V is conserved, equation ( |2.22j ) transforms to 



n 



[2.23) 



i=l 



where s = S/N^ and n = N-q/V stands for the entropy per baryon and the baryon number 
density. The comparison between equations (|2.12|) and ( |2.23|) gives the variation rate of 
entropy per baryon 



nTs = dii 1 (p {2) - pf ) 



i=i 



C+^-|2r 7 )^ + 2, / {x.x) 



(2.24) 
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3 Exact solutions 



To integrate equation ( [2.20D one needs a second set equations of state, involving the bulk 
viscosity coefficient C and the shear viscosity coefficient rj. Herein, we suppose 



K L) — 1 



—Toe \ 



2k 



(3.1) 



where Co > and rjo are constants. When the cosmological model is isotropic, i.e. x % = 
jo/(D — 1), i = l,...,n, and the shear viscosity is not significant, the expression H = 
7 e~ 7 /(Z} — 1) is the Hubble parameter. Then we get from equation (|3.1|) : C = ^^-( H, i.e. 
the bulk viscosity coefficient is linear proportional to the Hubble parameter. Such kind of 
the second equation of state describes the so called "linear dissipative regime" in the FRW 
world model (see, for instance, 0). As we study an anisotropic cosmological model, we 
must involve a shear viscosity as well as a bulk one. So, we propose equations ( |3.1| ) for 
the anisotropic model as a simplest generalization to the of the second equation of state 
describing the linear dissipative regime. 

In order to integrate equation ( |2.20|) , we use the orthogonal basis 



u 



(u, u) 



,f 2 , 



f n G JR n , 



(3.2) 



where the vector u was introduced by equation (|2.18|) The orthogonality property reads 

(u,fj) = 0, (^,4) = S jk , (j,k = 2,...,n). (3.3) 

Let us note that the basis vectors f2, . . . , f n are space-like, since they are orthogonal to the 
time-like vector u. The vector x e IR™ decomposes as follows 



u 



7o 



(3.4) 



Hence, under the above assumptions equation ( 2.20|) reads in the terms of coordinates in 
such basis as follows 



7o 



1 



V- Co) 7o -7 



7o 



D 



D-2 



k 2 A {h^ - h^) e ^ (1) - 2 ho+27 _ fc« E « =2 (x, ftf 
(x, f,-) + [(1 + ^0)70 - 7] (*,^> 
The integration of equation (|3.6|) gives 

(x,fj) =^e 7 - (1+,?oho , 

where p J is an arbitrary constant. 

Further we determine the time gauge by 



(j 



(3.5) 
,n). (3.6) 

(3.7) 



7 = fc 7o 



(3i 
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where k is a constant. 

By substituting the functions (x, fj) into equation ( |3.5|) we obtain the following integrable 
by quadrature ordinary differential equation 



70 - \k - 1 + 



+ Co 7o 2 



D-2 



,2(fc-l) 70 



D- 1 

In what follows we accept the agreement 



,J 1 g-2»7070 



3=2 



dt > 0, 



(3.9) 



(3.10) 



i.e. the cosmic time increases during the evolution. 

It should be noted that the solutions of equation (|3.9|) corresponding to different sets of 
the parameters h^ l \ h^ 2 \ Co and i] may lead to nonsatisfactory from the physical viewpoint 
cosmological evolutions. Further, we study only ones, which satisfy the following consistency 
condition: neither the density p^ 2 '(t) nor the variation rate of entropy s(t) have negative 
values on any time interval. 



3.1 The isotropic model 

The isotropic model is described by the metric given by equation fl2.1|) with 



a = e 



** = e WP-i) j = i n . 



(3.11) 



The scale factor a of the whole isotropically evolving space M\ x . . . x M n can be obtained 
by integration of equation (|3.9|) with 



= 0, j = 2, . . . ,n. 



(3.12) 



It can be proven that the above mentioned consistency condition leads to the following 
constraint 

h m _ h V) _ 2C > 0. (3.13) 

At first we present the special solution to equation ( |3.9| ). The special solution describes the 
asymptotical behaviour of all solutions at late time. It is the steady state solution 



a ~ exp 



2A(l + £o ) 



-t. 



\ (D — 1)(D — 2) 
for = 2 and shows a power-law behaviour 

t 2/[(2-fcW)(£>-l)] j (1) _ t 



« =A = A//t 2 



a ~ 



(3.14) 



(3.15) 



for 7^ 2, where t s is the synchronous time. Moreover, we present the deceleration 
parameter 

4S -2(*Zj*>=-l + £zi(2-fce>), (3.16) 



a 2 



the density ratio 

and the overall pressure 



P (2) /P (1) 



p (l) + p(2) 



^0 



2Co 



^)(l+5 )p (1) . 



(3.17) 



(3.18) 



The variation rate of entropy is positive (if Co > 0) and nTs ~ p^\ 

To obtain the general solution to equation (|3.9|) in the isotropic case we suppose that the 
parameter k specifying the time gauge by equation (|3.8|) is 



k = 1 - h?>/2. 



(3.19) 



(We note that the time t becomes synchronous if the perfect fluid component appears as 
A-term, i.e. = 2). This yields that equation ( |3.9| ) under the condition given by equation 
( 3.12|) is integrable by elementary methods. The result is 



i/[(D-i)(hW-hW-2<; )} 



a 

P (1) 
Q 

P (2) /P (1) 



a 



X 1-T 2 J 

Aa (n-i) [h m-2] 



■1 + 



D 



2 _ fc(D + _ h ® - 2Co) (1 - r 2 ) 



[l+e )/r 2 -l, 



:i + eo) 



1 - /> (2) - 2Co 



+ 1 - 



where we introduced the following variable 



r = tanh 



k 2 (L>- 1)A 
\2(D-2)(l + e ) 



(t - 1 ) 



t > t 



o- 



(3.20) 

(3.21) 
(3.22) 
(3.23) 
(3.24) 

(3.25) 



Now we analyze the role of viscosity in this isotropic 2-component model using the ob- 
tained exact solution. The main features of the model are the following. Under the weak 
energy condition (p^ + p^ > 0, p^ + p^ + + > 0), which leads to the following 
restriction on the parameters 2 > > + 2(o, the Universe expands eternally from the 
initial singularity. Near the singularity we obtain in the main order a ~ pVK 2- ' 12 - 2< ^o)(D-i)] 
There is only nonsingular solution given by equation ( |3.14|) . 

All solutions describes the accelerated expansion at least on the late phase of evolution 
if > 2(D - 2)/{D - 1). The solutions given by equation ( [3.20|) describe the period of 
decelerated expansion if /i (2) +2Co < 2(D-2)/(D-l). We note that the cosmic deceleration 
phase is important within mechanism of the clumping of matter into galaxies (see, for in- 
stance, @). Under the assumption 2 > > 2(D - 2)/(D - 1) > +2( the decelerated 
expansion takes place during the time interval (to,t*), where t* is defined by 



arcoshwl + 



t*-tn 



2(D-2)/(D-l)-fe( 2 )-2Co 



fe(i)-2(£)-2)/(D-l) 



2(D-2)/(D-l)-fe( 2 )-2Co 
h(l)-2(D-2)/(D-l) 



(3.26) 
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Equation ( |3.26|) shows that introducing of the bulk viscosity reduces the phase of the decel- 
erated expansion. 

The density ratio given by equations ( |3.17| ) and (|3.23|) exhibits the following property 



lim 



(2) /p (D 



2Co 



-/i(2) _2Co' 



(3.27) 



So, the bulk viscosity allow to resolve the coincidence problem which appears in the unviscous 
2- component model. 



3.2 The anisotropic model 

Now we consider the general anisotropic behaviour of the model. Once integrating equation 



( |3.9| ) we get the first integral of the form 7 = -F(7 ,C), where F is some function and 
C is an arbitrary constant. Substituting 7 and the functions (x, fj) given by equation 



^T7| ) into equations ( p.21| ) and ( |2.24| ) we express p^ and nTs via 70 . The subsequent 



analysis of the expressions for 70, p^ and nTs shows the presence of solutions with physically 
nonsatisfactory behaviour near the initial singularity. To exclude such solutions we put the 
following restrictions on the parameters 

2 > h (1) > h {2) + 2Co > -2r] > 2Co, -2r? > h {2) > 0. (3.28) 

These restrictions guarantee the following properties for all solutions: the density p^ and 
the variation rate of entropy s are positive during the evolution, which starts at the initial 
singularity of the Kasner type and proceeds eternally to the subsequent isotropic expansion; 
the density ratio p^/pW tends to the nonzero constant as t — ► +00 (see equation (|3.27| )). 



To study these behaviour in detail let us obtain the exact solution. We note that the 
equation 7 = F(jo, C) is integrable by quadrature for arbitrary parameters h^ 1 ', h^ 2 \ Co and 
i]q. To express the exact solution in elementary functions we put the following relation on 
the parameters 

h W _ 2h {2) - 2r] - 4Co = 0. (3.29) 

For the nonviscous model (r/o = Co = 0) the relation reads = 2h^ 2 \ The latter cor- 
responds, for instance, to the so called ACDM cosmological model with A-term [h^ = 2) 
and dust [h^ = 1). We remind that the parameters obey the inequalities given by formula 
p.28| ). Comparing equation ( |3.29| ) and formula ( |3.28| ), one gets > 0. 



Now we start the integration procedure in the time gauge defined by equations ( p.8|) and 
( |3.19| ) . The integration of equation ( |3.9j ) gives 



e*»M = C T(sina + T 2 COSa) , (3.30) 

1 — T 2 

where we cancelled the constant Sj=2(p') 2 by introducing the following constants 



C > 



2/i(2) " . \l ^A{hW-hW) Sj=2(P jN 



J2(P j ) 2 , oc = arcsin -2 G [0, vr/2], (3.31) 
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the parameter (3 is defined as follows 

P = h <X) - h &) - 2Co > 0. 



(3.32) 



The variable r was introduced by equation (|3.25 ). Substituting equation (|3.30| ) into equation 
(|3.7| ) and taking into account equations fl3.8|) and (|3.19|) , we obtain by integration 



(D - 2)/3 



-p 3 In 



r 



sin a + r cos a 



V/3 



(3.33) 



where g J are arbitrary constants. Substituting equations (J3.30|) and (|3.33|) into the decom- 
position given by equation Q3.4p, we get 



x = In 



Co 



v sin a + t cos a) 



V/3 



u 



(u, u) 



+ In 



r 



r(sina + r cos a) 



VP 



r + s, 



(3.34) 



where the vectors r G H n and s G H n are defined as follows 



E 

i=l 



r e, 



Owing to the orthogonality property given by equation the coordinates r l and s l of 
these vectors in the canonical basis e\, . . . ,e n satisfy the following constraints 



E^ = ( u > r ) = i, E rf * ( r = ( r > r > - E d A 



i=l 



i=l 



— (3.36) 



(g ~ g)g 
(D - l)hM D 

n 

dis 1 = (u, s) = 0. (3.37) 

i=i 

The constants r* may be called Kasner-like parameters, because of the existence of the 
constraints given by equation (|3.36j ). We remind, that the coordinates of the vector x in the 
canonical basis are the logarithms of the scale factors ai = exp[x*]. 
Finally, we present the exact solution: 



ds 2 



Co 



r(sin a + t cos a) 



1-r 2 



dt 2 + 



Cn 



(sin a + r cos a) 



2/(D-l)P 



X 



E 

i=i 



sin a + t cos a 



2r*/P 



e 2s 'ds 2 



(3.38) 



P 



T(sm a + t cos a) 



1 -(2-/i( 2 ))//3 



0" 



1-r 2 



(3.39) 



P i2) /p (1) = (l + e )F 2 (r) 



(3.40) 
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nTs 



I8k 2 A 



D - 1 
D - 2 



r + 110(3 f 2 (^ 



where 



*i(r) 



|(1 + t 2 ) sin a + r cos a; 
rfsina; + r cos a) 



1 



r sm a 



'1 + r 2 ) sin a; + 2r cos a 



(3.4i; 



(3.42) 



5 ^ 5 



The solution has the following integration constants: Co > 0, a G [0, 7r/2], to, r 1 , . • ■ 
s 1 , ...,s n . The constants r* and s l satisfy the constraints given by equations (|3.36|) and 
( |3.37| ). Then the number of free integration constants is 2n as required. The limit for 
a — > +0 of this exact solution is the isotropic solution obtained in section 3.2. 

Before we start the studying of the obtained exact solution near the initial singularity 
let us remind the mult i- dimensional generalization of the well-known Kasner solution || . It 
reads (for the synchronous time t s ) as follows 



ds 2 = -dt 2 s + J2At 2 s £l dsl 



(3.43) 



i=l 



Such a metric describes the evolution of a vacuum model under consideration. The Kasner 
parameters e % satisfy the constraints 



5> 

i=l 



i=i 



(3.44) 



The generalized Kasner solution describes the contraction of some spaces from the set 
Mi, . . . , M n and the expansion for the other ones. According to equation ( |3.44| ), the number 
of either contracting or expanding spaces depends on n (the total number of spaces) and 
di=\ >n (their dimensions). 

We note that the constraints given by equation (|3.36|) for r % coincide with these constraints 
for e' 1 when Co = "^o = and = 2h^ 2 \ i.e. the Kasner-like parameters r l become exactly 
Kasner parameters e % in the absence of viscosity. Therefore, if the parameters Co and 770 
are small enough, then the model describes a behaviour of the Kasner type as r — > +0, i.e. 
towards to the initial singularity. However, too strong viscosity suppress a behaviour of the 
Kasner type. It can be shown that if the parameters Co = ^0 = are large enough, then the 
model describes expansion of all factor spaces Mi, . . . , M n near the initial singularity. 

One can prove that the final stage of the evolution (r — > 1 — 0) exhibits the isotropic 
expansion. The asymptotical behaviour of the model for r — > 1 — is described by the exact 
solution given by equation ( |3.14j) - (|3~T8|) . 
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